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1 Introduction 

In this paper we solve the problem of classification of the consistent pairs of the lattices 

u x = F(ut,u,v), v x = G(u,v, u_i), (1) 
u y = P(u-!,u,v), v y = Q{u,v,v{). (2) 

Here u = u(n, x,y), v = v(n, x, y), n E Z and subscripts denote both the partial derivatives with 
respect to x, y and the shift with respect to n. The nondegeneracy conditions are assumed 

PvPuiG u G v _ 1 P v P u _ 1 Q u Q Vl ^ 0. (3) 

Due to the condition F v ^ 0, the first equation of the lattice (pQ) can be solved with respect to the 
variable v, then the second equation of the lattice rewrites as a lattice of the Ruijsenaars-Toda 

type una Gam 

Uxx = A( u l,x,Ux,U-l,x,Ul,U, U-l) (4) 

while its symmetry takes the form 

u y = B(u x ,ui,u,u-i). (5) 

Clearly, the roles of x and y may be interchanged and equation (j3]) may be replaced by a lattice 
of the form 

u yy = C(ui : y,u y ,u-i ty ,ui,u,u-i). (6) 

The classification problem for the lattices (|4|) was solved, in a different setting, in our paper [5], 
see also [U El [9] . The approach based on the relation to the lattices ([TJ allows to reproduce 
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this result in a slightly more general way and is promising for further generalizations. The 
classification of the integrable lattices of the form (HD was obtained by Yamilov [T3] . The main 
postulate in his work was that, instead of ([2]), a symmetry of a high enough order exists, and 
moreover, the lattice was assumed to be Hamiltonian. Our aim here is to obtain the answer 
under the minimal restrictions. However, in comparison with the result of Yamilov, our list 
contains only two more pairs. These pairs are not Hamiltonian, and they are reducible in some 
definite sense. We discuss these curious examples separately in Section [5j 

The interest to the lattices ([TJ is explained by their close relation to many other integrable 
models, see e.g. [10 1 PT2]. The higher symmetries of such lattices generate the hierarchies of 
the evolution systems of the nonlinear Schrodinger type. In this context equation ([2]) defines 
the first negative flow of the hierarchy and corresponds to some hyperbolic system, see Sec- 
tion 14.11 Although not all integrable two-component hyperbolic systems can be obtained in 
this way, this correspondence is a source of important examples. Also, the linear combinations 
of the flows ([T]), ([2]) constitute a class of equations containing Ablowitz-Ladik and Sklyanin 
lattices [131 ITS] . 

In Section [4.3l we demonstrate that the transform to the Ruijsenaars-Toda lattice ([4]) exhibits 
a hidden discrete symmetry of the lattice ([1]). Namely, it turns out that the lattices of the 
form (J3J) for the variables u and v coincide (after an appropriate substitution) , and this allows to 
expand the equations onto the square grid and brings to the discrete Toda lattices [HI [T71 [H] . 
It should be noted that conversely, the lattices of the form ([!]) or (j3J) can be obtained from the 
discrete Toda type equations under the continuous limits [19] . 

Section [7] contains the list of the consistent Hamiltonian pairs ([I]), ([2|) and the corresponding 
lists of the hyperbolic systems and Ruijsenaars-Toda lattices. 

Returning to the setting of the problem, notice that it is natural to consider equivalent lattices 
related by the point substitutions 

u = 4>(u), v = ip(v) (7) 

the scalings x = ax, y = (3y and the renamings 

u «-> v, n «-> —n, F «-> G, P <-> Q, (8) 
x <-> y, n <-> — n, F P, G <-> Q. (9) 

We use these transforms in order to bring a lattice to a simpler form. Moreover, the substitution 

u = u + ax + j3y + 771, v = v + ax + (3y + 771 (10) 

is often useful. It is clear that it brings, in general, to a nonautonomous lattice. However, if the 
right hand sides of the lattices contain only the differences u — v, u±\—v then this transformation 
preserves the class under consideration. For brevity, we refer to all mentioned transformations 
as to admissible substitutions. Our main result is the following theorem. 

Theorem 1. The consistent lattices ([1]), ([2]), such that the nondegeneracy conditions ([3]) are 
fulfilled, are brought by the admissible substitutions to one of the Hamiltonian pairs in the list \7.1[ 

u x = a(u,v)5 v H, v x = -a(u,v)S u H, H = K (u n+ i, v n ) + L(u n , v n ), (11) 
u y = a(u,v)5 v R, v x = -a(u,v)5 u R, R = M(u n , v n+1 ) + N(u n , v n ) (12) 

(5 U H = d u ^2, n T n {H) denotes the lattice variational derivative), or to one of the pairs 



v + (1 + e)u + u\ 



v-i + {l + e'^v + u 
v-i — e~ 1 u 



ev — u\ 
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v + (l + e 1 )u + u_i vi + (1 + e)v + u . 

% = zr , = , (*7) 

e l v — v\ — eu 

u x = f(v-u 1 ), v x = ef(v-i - u), (X 8 ) 

u y =p(eu-i-v), v y =p(eu-vi) (Y 8 ) 

or to the linear lattices. 

2 Necessary conditions 

In this Section we deduce some consequences from the compatibility condition of the lattices 
under consideration. Let T denote the shift operator n — > n + 1. 

Proposition 1. If the lattices (|T]), ([2]) commute then a constant e and functions a(u,v), 
b(u,ui), c(v,v\), h(ui,v), r(u,vi) exist and are unique up to the scaling (a,b,c,h,r) — ► 
(ka,kb,kc,h/k,r/k), such that 

F U1 = ah, F v = bh, (13) 

G v _ 1 = -eaT~\h), G u = -T" 1 ^), (14) 

P u _, = -eaT-^r), P« = -T _1 (6r), (15) 

= ar, Q„ = cr. (16) 

Proof. Computing u xy in two ways gives the equation 

F U1 T(P) + F U P + F„Q = Pu^T'^F) + P U F + P V G. (17) 

Differentiating of this expression with respect to u\ and U-\ yields F UUl P u _ x = -P uu _ 1 i ? Ul , which 
implies 

F Ul = a(u,v)h(ui,v), P u _ 1 = a(u,v)r(u-x,v). 
Using the symmetry (|8|) gives analogously 

G v _ x = a(u,v)h(u, v-i), Q Vl = a(u, v)r(u, vi). 

Let us prove that, moreover, the following relations hold 

T(P V ) = -br, aF v = abh, (18) 

T~ 1 {Qu) = -cr, aG u = ach, (19) 

T"- 1 ^,) = -bh, aP v = abr, (20) 

T(G U ) = -ch, aQ u = acr (21) 

with some functions 

b = b{u,u\), c = c(v-i,v), 6 = 6(n_i,u), c = c{v,v\). 

In order to obtain the first line, differentiate (I17p with respect to v\. This yields F v Q Vl + 
F Ul T(P v ) = 0, or aF v /(ah) = —T(P v )/r. Here the left and right hand sides depend respectively 
on u\, u, v and u, u\, V\, so that their common value is some function b{u, U\). In order to prove 
the rest, it is sufficient to use the symmetries flHJ), ([9]). 
Now, one easily obtains from (I18h the relation 



a{u,v) b(u,u\) _ T{a(u,v)) c(v,v\) 



Ku,v) T(b(u-i,u)) T{a{u,v)) T(c(v- lt v)) 
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This implies a/a = tp(v)/cf)(u). One may assume a = a without loss of generality, taking into 
account the substitution 

a — ► aip, a — > a(f>, h^h/ip, r^-r/ip, h^h/cj), r—*r/(j). 

Then b = eT(b), c = eT(c) with some constant e and functions h, h and r, f are related by 
equations 

h = -eT- 1 (h), f = -eT~ l {r). 

Taking all together proves the required equations. ■ 

The system of equations (|13p - (|16p is overdetermined and further information about the form 
of the lattices will be obtained mainly by analysis of the compatibility conditions of this system. 
However, this cannot solve the problem completely, since the functions F, P and G, Q are found 
from this system at least up to the addition of arbitrary functions on u and v respectively. 
Therefore, after some point we will need additional information. It turns out that use of the 
following proposition allows to determine the right hand sides of the lattices up to few arbitrary 
constants. The final answer is obtained then by the intermediate check of the consistency. 

Proposition 2. If the lattices ([I]), ([2]) commute then the following equations hold 

h v Q + h Ul T(P) = -h{T(P u ) + Q v + k 1 ), a u P + a v Q = a{P u + Q v + k 1 ), (22) 
r u F + r Vl T(G) = -r(F u + T(G v ) + k 2 ), a u F + a v G = a(F u + G v + fc 2 ), (23) 

where k\, k 2 are some constants and the functions a, r, h are defined in the Proposition^ 
Proof. Equations (|22p . (|23p are equivalent to the conservation laws 

D y {\ogF ul ) = (1 - T)(P U ), D y (logG v ^) = (1 - T^XQ,), (24) 
A^logiVj = (1 - T- 1 )^), D y (logQ Vl ) = (1-T)(G V ), (25) 

where the left hand sides are replaced in accordance to the formulae (|13|) — (| 16 [> . In its turn, the 
first of these conservation laws is obtained by differentiating of (fTT|) with respect to u\, and in 
order to obtain the other ones it is sufficient to apply the symmetries (JSJ) , (EJ) - B 

Remark 1. We see, comparing the equations and (I13p . (|14p that for the Hamiltonian 
lattices h = H U1V and e = 1. In the general case parameter e remains undetermined almost 
till the end of calculations. This makes necessary to consider several additional branches with 
£ / 1, leading to essential complication of the analysis (see Propositions [H below). Only in 
two of these cases the solution is not empty. Moreover, it is not difficult to prove that in the 
Hamiltonian case the second column of the equations (|22p . (|23p turns into identities and the 
constants fei, k 2 vanish. 

In what follows we use the substitutions ([7]). They act on the functions under consideration 
in accordance to the rules 

F = <f>'(u)F, G = i>'{v)G, P = <j)'(u)P, Q = ^'(v)Q, 

h = VW>WM), r = r/(0'(«y («i)), 

a = <fr' (u)ip' (v)a, b = eft' (u)(J)' (ui)b, c = ip' (v)ijj' (vi)c. 
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3 Proof of the classification theorem 

We will see soon that it is convenient to divide all lattices into two classes, depending on whether 
a is of the form a(u, v) = a\{u)a2{v) or not. We start from the case when this factorization does 
not hold. The form of the lattices from this subclass (which contain the larger part of the list) 
is defined more exactly in the following proposition. 

Proposition 3. If the lattice dU) satisfies the equations ()13p . (|14p at aa uv ^ a u a v then it can 
be brought, by a substitution ([7]), to one of the following types: 



\u x = aui + f\(u)v + f(u), 
\v x = -eav-i + gi(v)u + g(v 
( a a. 



i 



ayitV, (26) 



i,j=0 



Ux = ^ + /(«), 2 

V — U\ I 

£CL &Oj U 



Oij-uV. (27) 

i,j=0 



v-i — u 2 

Proof. Cross-differentiating of the relations (|13j) yields 

a v h + ah v = b Ul h + bh Ul . (28) 
Dividing by h and applying the operator d v d Ul — (logh) VUl yields 

aN Ul = bM,„ N = ^ - S|, M — KlUl 3kui 



h 2/i 2 ' h 2h 2 



Differentiating with respect to u gives a u N Ul = b u M v . Since aa uv / a u a v , hence a u b ^ ab u 
and therefore M v = N Ul = 0, that is M = M(u±), N = N(v). It is easy to obtain the rule of 
transformation of these functions under the substitution ([7]): 

W ' ^ ft 2(0') 2 ^ 2(^') 2 

We set M = N = by an appropriate choice of 4>, tp. Denote h = H~ 2 , then these equalities 
and equation (|28l) take the form 



H vv — 0, H UlUl — 0, a v H — 2aH v — b Ul H — 2bH Ul . (29) 

Notice that we may still apply the Mobius transformations of v and u\ . This allows to bring H 
to one of the forms H = 1 or H = v — u\ (depending on whether the original polynomial is 
reducible or not). 

If h = 1 then the system (fT3"|) takes the form F Ul = a(u,v), F v = b{u,u\), whence 

F = a{v)vu\ + (3(u)v + ^(ujui + S(u). 

Analogously, 

G = x{v)v-\u + X(v)v^\ + fi(v)u + v{v), 

moreover — ea = G v _ 1 = —eF Ul , and equations (f2"Uj) follow. 

If h = (v — tii) -2 then one easily finds from equations (|29p that a = a2V 2 + a±v + ao, 
6 = — a2u\ — a\U\ — ao, where the coefficients depend on u. Solving the system (fTBI) yields 

a(ii)vMi + P(u)v + 7(14)1*1 + 



v — ui 
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where a = a>2, j3 + 7 = a%, 5 = ao- Analogously, 

_ x{v)v-\U + \{v)V-\ + fi(v)u + v{v) 

C / = 

u_i — u 

and one proves that a is a biquadratic polynomial by comparing the formulae 

a = a(u)v 2 + (P(u) + j(u))v + 5(u), ea = >c{v)u 2 + (A(«) + fj,(v))u + u(v). 

The lattice takes the form (j27|) with 2/ = (3 — 7, 2g = A — fi. ■ 

The similar statement is valid for the lattice ([2]) as well, due to the symmetry ([9|). However, 
the substitutions © for both lattices may be different, so that the second one should be written 
for some another variables U = U(u), V = V(y): 



U y = -eAU^+ Pl (U)V +p(U) 
V y =AV 1 +q 1 (V)U + q(V), 

£ A eAy / t t \ 



1 

A = AtjU'Vi, (30) 

i,j=0 



A Au 
U-Vi 2 



V y = --£- + q(V), iJ=0 



A=Y J A ij U i V^. (31) 



The lattices (|2T)j) . (|27|) and ([5U|) . ([3T|) combine into three variants of the pairs, up to the symmet- 
ry 0. In each case some additional information about the right hand sides can be obtained by 
comparing the formulae for the functions a, b, c corresponding to both lattices. The following 
lemma is useful as well. Here we denote 

F(U 1 ,U,V) = U'F(u 1 ,u,v), G(U,V,V. 1 ) = V'G(u,v,v. 1 ) 

for the right hand sides of the lattice ([T]) rewritten in the variables U, V. 

Lemma 1. // the lattice ([T]) is of the form (I26j) then its symmetry (|2|) satisfies the conditions 
Puu = Qw = 0. If the lattice JT]) is of the form <\27\) then P uuu = Q vvv = and, moreover, 
P and Q coincide as the functions on their arguments, that is P(u,v,v\) = Q(u,v,v\). 

Analogously, if the lattice ([2|) is of the form (j30|) then Fuu = Gyy = 0, and if the lattice (|2|) 
is of the form flE) then F uuv = Gyyy = and F{U U U, V) = G(U X , U, V). 

Proof. Let us use the first of the relations (|22p . If h = 1 then the equalities P uu = Q vv = 
follow immediately. If h = (v — u\)~ 2 then differentiating with respect to u\ gives 

Q = T{P) + (v- ux)T(P u ) + i(« - Ul ) 2 T(P uu ). 

Differentiating once more yields P uuu = 0. Moreover, the last equation is nothing but the Taylor 
expansion in the second argument for the function P(u, v,vi). The equations for F, G follow in 
virtue of the symmetry ([9]). ■ 

Proposition 4. The consistent lattices with the condition aa uv 7^ a u a v are brought by the 
admissible substitutions to one of the pairs (X± , Yj )-(^4 . ) or ( IX7HY7D . 

Proof. 1) Let the lattices be of the form (|27p . (|3ip . Then one can set h = (v — u±)~ 2 and 
r = (U — Vi)~ 2 U'V[, after scaling of y. The comparison of the functions a, b, c gives 



MU,V) , , . AjU^U) A(V U V) 
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where 

2 2 

a(u,v) = o<juV, A(U,V) = Y AijTPV 3 . 

i,j=0 i,j=0 

The differentiation of the first equation yields 

{a u U' + aU")V = A V U', (a uu U + 2a u U" + aU'")V = Auu(U') 2 + A V U" , 
(3a uu U" + 3a u U"' + aU IV )V = 3AuuU'U" + AuU'" 

and elimination of Ajj, Ajju brings to the relation 

{2U"'U' - 3{U") 2 )U'a u = (4U"'U"U' - U IV {U') 2 - 3(U"f)a. 

Since aa uv ^ a u a v this implies 2U"'U' — 3{U") 2 = 0, that is U is a linear-fractional function. 
The same is true for V as well. 

Since the lattice (|3ip is form invariant under the substitution U = M(U), V = M(V) where 
M is an arbitrary Mobius transform, hence one may set U = u without loss of generality. Then 
one may apply the Mobius transformations to the variables: u = M(u), v = M(v). Under this 
transformation the function V is changed in accordance to the formula V = MVM -1 , so that 
it can be brought to one of the forms V{v) = 5v or V(v) = v + 25. The system (|32p becomes 
equivalent to the equation a(u, v) = eA(v, u) which defines A and the relation 

ea(u,v)V'(v) = a(V(v),u) 

which means that function a possesses some generalized symmetry property. All biquadratic 
polynomials which satisfy this identity can be found directly. They are listed in the Table [U (up 
to the inversion (u,v) — > (1/u, 1/v); the restrictions on the parameters are introduced in order 
to avoid the intersections and to provide the condition aa uv ^ a u a v ). Now we only have to find 
the coefficients of /, g, p, q which are quadratic polynomials, in accordance to the Lemma [TJ 
This is done by the direct computations which show that the most cases are empty. As a result, 
we obtain the pairs ( Xi|Yi ), ( X7IY7 ) and also the pairs 

a a 
u x = h u — v + (3, v x = h u - v + p, 

V — U\ V-l — u 

a a 

u y = — -u + v + j, v y = — - u + v + 7, (33) 

■u_i — v — 2d u — v\ — 2d 

where a = (u — v — 5) 2 + a, and 

a r n a in . \ 

u x = dv + pu, v x = h u + ip + a)v, 

V — U\ V-l — u 

8a r , . a 

u v = — + dv + (7 + a)u, v v = u + jv, (34) 

U-i — dv u — dvi 

where a = u 2 + auv + dv 2 . However, these lattices can be brought to the particular cases 
of the lattice flA7|Y"iD by the transformation (fTO]) . For the lattices ([33]) one should apply the 
substitution 

u — > u + (P + S)x + (7 — d)y — 6(n — 1), v — > v + (f3 + d)x + (7 — d)y — dn 



and the new variables satisfy the lattices ( X\ \Yi ) with a = (u — v) 2 + a. For (|34|) . the analogous 
substitution is conjugated by exponentiation. 
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Table 1. Solutions of the equation ea(u,v)V'(v) = a(V(v),u). 



e 


V(v) 


a(u, v) 




1 


V 


CL22U 2 V 2 + ai2Uv(u + v) + O02 {u 2 + V 2 ) 








+a\\uv + aoi(u + v ) + aoo 


(.Xi.Yil 


-1 


V 


(u — v)(auv + (5{u + v) + 7) 




1 


v + 25 


(u - v - 5) 2 + a, 5^0 




-1 


v + 25 


u — v — 8, 5^0 


— 


#1 


ev 


au 2 v 2 + eu + v, a/0 






—v 


auv(u + ev) + (3{eu + v), af3 / 




-l 


—v 


au 2 v 2 + (3{u 2 + u 2 ) +7, 07 / a 2 




l 


Sv 


u 2 + auv + (5u 2 , 5/0,1 


(IHI) 


-l 


Sv 


u 2 -^ 2 , 5/0,1 






e~ 2 v 


eu + v, e 2 / 1 


(X 7l Y 7 ) 



2) Let the lattices be of the form (ETJ) , §M) , h = (v - m)' 2 , r = U'V{. Then 
4 , A U UV + A 10 U + A 01 V + A 00 



i,i=o 

^nC/C/x+^oiCZ-piCf/i; 



&(it, ui) = —a(u, U\) 



U'U[ 



, > , x An^ + AioFi + giC^) 
c(u,ui) = -ea(u,ui) = . 

As in the previous case, the first equation implies that functions U, V are linear-fractional. The 
other two equations are equivalent to relations 

AnV^+Aw _ A u U(v)+A 01 

V'(v) ~ £ U'(v) ' [6b) 

Pi{U{u)) = A 01 V(u) + Aqq qi(V(v)) _ A 10 U{v) + A 00 

U'(u) V'{u) ' V'(v) 6 U'(v) ' 1 ' 

We assume, without loss of generality, that either A = UV + 1 or A = ell — V . Equation (|35p 
implies U(u) = cV{u)~ e (and, consequently, e = ±1) in the first case, or U(u) = V(u) + c in 
the second one. The appropriate substitutions, the Mobius one in the lattice (|2T[) and the linear 
one in (j30j) . reduce the problem to the following cases: 



e = l,U = u, V = v~\ a = -v(u + 5v), A = UV + 6, 5^0; 
e = — 1, U = u, V = v, a = A = uv + 1; 
U = u, V = v, a = A = eu — v + 6. 

Moreover, the functions p±, q\ are defined from the relations (|36p and functions f(u), g(v), p(u), 
qiV) from the equations 

/" = 0, ((VgY/Vj = 0, p"> = 0, {q/V') vv = 0, 

in accordance to the Lemma [TJ In the first case, after defining the coefficients by the direct 
computation and the substitutions u — ► e u , v — > e v and (fTO)) we obtain the lattices flA^IY^D . The 
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second case turns out to be empty and the third one contains, at e = 1, the pair ( X2II2 )> U P to 
the admissible transformations. 

3) Let the lattices be of the form ()26|) . (|30j) . h = 1 and r = U'V{. The comparison of the 
expressions for the function a brings to the equation 

(a u uv + ai u + a iv + a 00 )U'V = A U UV + A 10 U + A m V + A 00 . 

Computing aa uv — a u a v one obtains the relation 

(a u a o - ai a i)U'V' = AuA o - Ai A i 

which implies U" = V" = 0. Therefore, we may assume, after the linear substitutions, U = u, 
V = v, A = a. Lemma[T]says that the right hand sides of both lattices are linear in any variable. 
Several equations for the coefficients follow from the relations for the functions b, c: 

6(n, u\) = a\\uui + ao\U\ + fi(u) = eauuui + eaoiu — pi(ui), 
c(v, vi) = eauvvi + eai v - gi(vi) = auvv\ + a w v\ + gi(^). 

After the final check of the compatibility condition we arrive to the pair ( X4IY4 ). ■ 

In order to finish the classification we have to consider the lattices with the function a of the 
form a(u,v) = ai(u)a,2(v). 

Proposition 5. The consistent lattices with the condition brought by the 

admissible substitutions to one of the pairs ( LYstjTsI ), (JXfsJEI or ( Ell ED - 

Proof. The substitution ([7]) allows to reduce the problem (temporarily) to the case a = 1. 
Then the compatibility conditions for the systems (|13p -(|16 p are 

hv , hui hv r vi t ,T U fu r vi /o^7^ 

-7- = o«i + 0-7-, e—- = c v + c—, e — = b u + b — , — = c vi +c — . (37) 
h h h h r r r r 

As a corollary, the equations hold 

6«(log/i) OTl = 0, (log h) VUl = 0, b Ul (logr) UVl = 0, c v (logr) UVl = 0. 

We will use also the following consequences of the relations (f22j) . ([23)) : 

P u + Qv +h = 0, F u + G v + k 2 = 0. (38) 

1) At first, assume that both quantities (logh) vul and (logr) m)1 do not vanish. Then b, c are 
constants and be = s. A scaling of u and v allows to set b = —1, c = — e. One obtains, after the 
integration of the equations (fT3|) -([T6 |) and taking (f38|) into account, the lattices of the form 

Ux = f(v - ui) + aiu + a, v x = ef(v-i - u) + fcv, 
u y = p(eu^i - v) + 71U + 7, v y = p(eu - vi) + Siv. 

The direct computation proves (notice that f'p" 7^ by assumption) that the linear terms are 
zero and brings to the pair ( LYg|YgP . 

2) Now let either (logh) VUl = or (logr) UVl = 0. Taking the symmetry ([9]) into account, we 
assume for definiteness that h = m{ui)n[v). Then the two first equations (|37p take the form 

n' m' m! n' 

— = b Ul +b — , e — = c v + c—, 
n m m n 
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whence 

m = fim, n = vn, b Ul = v — fib, c v = e/x — vc. 

Here we have to consider several subcases. 

2.1) Let jiv ^ 0. After scaling, assume fi = 1, v = —1, h = e Ul ~ v . Then 

b = e~ Ul m(u) — 1, c = e v n(v\) — e, 
and two last equations (|37|) are reduced to relations 

m'(u)r + fh(u)r u = 0, n'{v\)r + n(vi)r Vl =0, r u + er Vl = 0. 
If m = n = then we obtain the lattice of the form dAg|Yg[ ). Otherwise, 

b = p e - Seu - Ul - 1, c = -/e v+Svi - e, r = Xe 8{£U - Vl \ 
where at least one of the coefficients /3, 7 is not zero. The solution of equations (|13p . (|14p is 

F = e" 1 " 1 ' - f3e~ Seu - v + f(u), G = ee u ~ v -' - je u+Sv + g{v). 

The functions P, Q are easily found as well. One easily proves, by use of (|38p. that 5 = — 1, 
e = 1, f3 = 7 and the functions f,g are linear. After the direct computation of the constants 
and the substitution e u — > u, e - '" — > w (of course, it spoils the gauge a = 1) we come to the pair 

2.2) Let fi = 0, f 7^ 0. Taking 1/ = 1, /i = e" we obtain 

b = u\ + fh{u), c = e~ v n(vi), 

r u = 0, n'(vi)r + 71(1)1)^ = 0, er Vl = fh!(u)r. 

From here, b = u\ + eJu + /?, c = 'ye - '" - ' 5 " 1 , r = Xe Svi and the solution of equations (p~3|) . (fI3]l is 

F = (m + 5eu + /3)e" + /(u), G = -ee v ~ l - jue~ Sv + g(v). 

The second equation (|38p now reads fee" + /'(it) + 7^e" 51 ' + g» + k 2 = 0. Since 7 7^ in 
virtue of the nondegeneracy condition Q, hence 5 = 0. Taking r = 1 and solving equations 
([T5]). (fTBj) one obtains 

-P = —eu-i — (u + 0)v + p(u), Q = vi + 7ue - " + q(v). 

However, then the first equation (|38p implies 7 = 0. Therefore, this case is empty (as well as 
the case fi / 0, v = 0, due to the symmetry ©). 

2.3) Finally, let /x = v = 0. Then b = b(u), c = c(v±) and one can take h = 1. Then 
F = u\ + 6(n)v + /(u), G = — — c(v)u — g(v) and substitution into (f38|) yields 

6'(u)i; + /'(u) - c»it - + k 2 = 0. 

From here, 

6(it) = cm 2 + 2f3\u + /3, c(t>) = aw 2 + 271 1> + 7. 
Moreover, elimination of r from the equations (|37j) brings to the relation 
eft" + (6 /2 - bb")c = e 2 c" + e{c' 2 - cc")b, b = b(u), c = c{ Vl ). 
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Denote A& = b' 2 — 2bb" = 4(/3^ — a/3), and analogously for c, then this relation takes the form 

2ea + (A b + 2ab)c = 2e 2 a + e(A c + 2ac)b. 

If a 7^ then it follows that e = 1 and A b = A c = 0, that is b and c are full squares. The 
linear change of variables brings them to the form b = u 2 , c = vf. Then one finds from (|37|) 
that r = (uv\ — 1)~ 2 and solution of the equations (fT5j) . (fT6|) is P = U-\/(u-iv — 1) +p(it), 
Q = —vi/(uvi — 1) + q{v). Equations (f38l) say that functions f,g,p,q are linear and after the 
direct computation of the constants we obtain the pair ( A^ll^D - 



If a = then = 0, $7 = £7^. Since functions b, c do not vanish in virtue of ([3]), hence 
/3i = 7x = 0, that is the lattice ([I]) is linear. Direct calculations prove that either the second 
lattice is linear as well, or we arrive to a particular case of the pair ( bYg|Yg[ ). ■ 



4 Associated equations 

4.1 Hyperbolic PDE systems 

In accordance to the nondegeneracy conditions the equations (pQ), ([2]) can be solved with respect 
to the variables u±i, v±i: 

Ul = F(u x ,u,v), v-i = G(u,v,v x ), 
u-i = P(u y ,u,v), vi = Q(u,v,v y ). 

Therefore, these variables can be eliminated from the expressions for the mixed derivatives and 
some system of partial differential equations appears. The original lattices define its Backlund 
auto-transformation. The form of the system is given by equations 

Uxy = fiUxUy + fsU x + /2% + fi, v xy = g^v x v y + g 3 v x + g 2 v y + 01, 

where the coefficients depend on u, v. Indeed, consider the equalities 

Uxy = F Ul T(P) + F u u y + F v v y = P u _ L T~ l (F) + P u u x + P v v x . 

In the first one the variables u%, v\ should be eliminated only and we see that the expression 
for u xy does not depend on v x and is linear in u y . Analogously, elimination of U-i, v-\ in the 
second equality proves that the expression for u xy does not depend on v y and is linear in u x . 
The formula for v xy is proved analogously. 

Proposition 6. In virtue of the lattices from the list \7.1\ the variables u, v satisfy the respective 
systems from the list 17.21 

4.2 Ruijsenaars Toda lattices 

As we have already explained in the Introduction, elimination of the variable v allows to rewrite 
the lattice ([I]) in the form of Ruijsenaars-Toda type lattice and its symmetry ([2]) in the 
form ([5]). Moreover, the equation (|6|) is fulfilled as well. The triples of the lattices corresponding 
to the list 17.11 are enumerated in the list 17.31 The remarkable property is that elimination of u 
instead of v brings to the same equations. This can be proved without calculations since it is 
easy to see that the lattices from the list 17.11 are invariant with respect to the involution 

u n <-> av- n , x^-x, y -> -y, 
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and the lattices from the list 17.31 are invariant with respect to the involution 



-x. 



y 



-y, 



where a = 1 for flXiD , flYiD , Q-RiD and a = — 1 in other cases. It should be noted that this 
property is not invariant under the general substitutions ([7]). For an arbitrary choice of the 
variables u and v the corresponding Ruijsenaars-Toda lattices coincide after a suitable point 
substitution. As to the lattices QA^D , ( [T^j ), flAg| ), ( |Tg| ), they possess the above symmetry only if 
s= 1. 



4.3 Discrete Toda lattices 

Coincidence of the lattices (|4|) for the variables u and t> allows to expand the equations on 
the square grid and leads to the discrete Toda lattices. More precisely, this property means 
that the equations ([I]), ([2]) define the Backhand auto-transformations for the lattices (J1J), © 
correspondingly. This allows to introduce the variable u(n, m, x, y), so that the pair (u(n),v(n)) 
is identified with {u(n,m),u(n,m + 1)) at arbitrary m. Let us rewrite the lattices from the 
list 17.11 denoting the shifts by the double subscripts: 

u x = F(n ,i, u, ui >0 ) = G(u-x t0 ,u, u _i), (39) 
u y = P(u 0) -i,u, ui )0 ) = Q(u- 1>0 , u, u ,i). (40) 

Then the following properties are valid: 

1) equations F = G and P = Q are equivalent and can be rewritten as the discrete Toda 
lattice 

f(u ,i,u) + f(uo-i, u) = g(u 1)0 ,u) + g(u-ifl, u); (41) 

2) equation (|4ip is consistent with the dynamics on x and y, that is, the equations obtained 
by differentiating of this equation in virtue of ([39|) or ([4"0j) are its consequences; 

3) the variables u(n) = u(m,n) for any m satisfy the Ruijsenaars-Toda lattices (|4"|), 

4) the variables u{m) = u(m, n) for any n satisfy the analogous lattices (possibly, with the 
different right hand sides); 

5) the variables along the diagonals u(n) = u(m — n, n) satisfy a Toda type lattice with 
respect to the ^-derivatives 

Uxx = A( U x,Ul,U,U-l), 

and the variables along the complementary diagonals u(n) = u(m + n, n) satisfy a Toda type 
lattice with respect to the y-derivatives 

u yy = C(uy,ui,u,u-i). 

6) all these equations are Lagrangian. 



5 Non-Hamiltonian lattices 

The pairs ( X?\ Y7D and dATgjlg ) are of interest as counter-examples to the statement that existence 
of one symmetry implies the existence of the whole hierarchy. These lattices do not possess the 
higher symmetries at arbitrary /, g and e. In particular, the necessary integrability conditions 
fail for the associated lattices of the form In accordance to [13] the simplest of these 

conditions reads 



D x {logA Ulx ) elm(T-l). 
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The direct computation proves that elimination of v from equations (Xj) brings to the lattice 
(u x + l)(eu x - 1) / ui >x + 1 eu_i, x - 1" 



1 + e \u\ + eu u + eu^i 

and the above condition is fulfilled for this lattice only if e = 1. 

It turns out, however, that the nature of these counter-examples is not too profound. It is 
easy to see that the non-invertible substitution U = v — u\, V = eu — v\ reduces the lattices 
( LYg|YgD to the disjoint pairs 

U x = ef(U-i) - f(Ui), U y = 0; V x = 0, V y = ep(F_i) - p{V x ). 

Analogous, but not so evident substitution exists for the pair ( IX7II7P as well: the variables 

u= WM , v= d + e )(u + ) 

(£V-i — u)(ev — ux) (eii-i — v)(eu — v\) 

satisfy the equations 

U x = U(U 1 + U-e(U + U- 1 )), U y = 0; 

V x = 0, V y = -V(Vx + V - e(V + VLi)). (43) 

In both cases the quantities U (n) are the local first integrals for one lattice in the pair, and V(n) 
for another one. In these variables the consistency of x- and y-flows become trivial and is not 
related to integrability (although at e = 1 and for some /, p equations may "accidentally" 
coincide with integrable Volterra type lattices |23j). 

The fact that the local first integrals satisfy the closed lattice with respect to the second 
independent variable can be easily explained in the framework of the associated hyperbolic 
system. For the pair QX7II7P it is of the form 

_ K + 1)K + X ) „ _ (Vx-l)(Vy-l) 

Uxy ~ u + v ' xy ~ u + v 1 ' 

(notice that the units may be dropped, due to the substitution u ^ u — x — y, v ^ v + 
x + y, however this leads to nonautonomous lattices). This is the Liouville equation type 
system [20^ |2T| [22] . Elimination of the shifts in the first integrals (142 p leads to the invariants 
of this system, that is, to the quantities I(u, v, u x , v X} ■ ■ ■ ), J(u, v, u y ,v y , . . .) which satisfy the 
properties I y = 0, J x = in virtue of (j44]). It can be proved that, in the case of the hyperbolic 
system of the second order, all solutions of the equation I y = are expressed through the x- 
derivatives of at most two basic invariants. Therefore, the invariants U±, U, U—\ must be related 
by some differential constraint, as the first equation ([5]) shows. 
Finally, notice that the system (|44h implies 



(log U xy ) X y = 0, (lOg V xy ) x y = 0. 

This gives the formula for the general solution 

(U x + l)(Uy + 1) 



u = a(x)b(y) + c(x) + d(y), v = -u + 

under assumption a'b' ^ 0, and this ansatz reduces the lattices ( X7IY7 ) to the recurrent relations 

ec'-l d' + l 

a\ = ; — , c\ = aa\ — ec; Oi = — - — , d\ = e{bb\ — d). 

a' b' 

These relations are equivalent to some Volterra type lattice, for example the substitution a — > e n a 
and scaling of x brings to the equation (known to be integrable at e = 1 |23j) 

- 1 

£ 0"n,x 

O-n+l — O-n-1 
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6 Concluding remarks 

The integrable lattices of the form jl]) or admit the very important generalization: their right 
hand sides may contain parameters which depend on n and define the discrete spectra added by 
the iterated Backlund transformations. Such generalizations are known, actually, for all lattices 
from the above lists, see e.g. [121 E]- It turns out, however, that n-dependent lattices (pQ) are not 
consistent with the symmetries of the form ([2]), only with the higher symmetries of NLS type. 
The nature of this phenomena is not well understood for the present. It takes place in another 
situations as well, for example, the dressing chain v n>x + v n +i jX = (v n — i> n +i) 2 + a n is consistent 
with equation v n>y = (i> n +i — i^-i)" 1 only if a n = a n +i [21]. However, the exceptional pairs 
( X7IY7 ), ( Xg|Yg ) may depend on n, for example, the following pair is consistent: 



v n + (1 + e n )u n + u n+1 v n -i + (1 + e n \)v n + u r 



£ n v n — Un+i v n _i - e n \u r , 



_ Vn + (1 +£ n - 1 )u n + U n -i _ V n+l + (1 + e n )v n + U n 

^n,y — _i > Vn,y — 

£ n -l V n - Un-1 V n+1 ~ ^nUn 

Another sort of the nonautonomous generalizations is related to the master-symmetries of 
the lattices (fTjl. ([2]). see [T41I25] , 

The multifield lattices also should be mentioned. Probably, the earliest example of such kind 
was found in (2S1 EZ]. It generalizes the pair ( bYglYeD hi the rational form: 



u x = u\ + (it, v)u, —v x = V-i + (it, v)v; u. 



(u-i,v)-l y (u,vi)-l 

The vectors it, v satisfy, in virtue of these equations, the system 

U X y = {U,v)lly + ({Uy,v) ~ l)u, V X y = ~{u,v)Vy ~ ({U,Vy) + l)v . 



However, an analog of the lattice (Rq) is absent in this example. 



7 Lists 

7.1 Consistent pairs of the Hamiltonian lattices 

a a v a a u 
u x = -, v x = h— , (Xi) 

V — U\ I V—l — u 1 

a <i v a a u 

u y = + tt, v v = ir, \Xi) 

it_l — v 2 it — v\ 2 

a = a\u 2 v 2 + a.2iiv(u + v) + a^,{u 2 + v 2 ) + a^uv + 0:5(1* + v) + ate, 
H = log(iti - v) - - log a, R = - log(u - v x ) + - log a, 

u x = (it - ui)(u - v), v x = (v- 1 -v)(u-v), (X 2 ) 

It — U-l V\ — V 



v — lt_i m — u 

a = u — v, H = uv — uiv, R = log(it — vi) — log(it — v), 



2) 



Ux = {l + e u ^ u ){l + e u ~ v ), Vx = (l + e v ~ v -i)(l + e u - v ), (X 3 ) 
1 + e u_ 1 -u l + e v ~ Vl 

U y ~ 1 _ e u.!-v ' V V ~ l _ e u- Vl ' ( Y 3> 

a = l + e v ~ u , H = -e Ul " v - e u ~\ R = log(l - e vi ~ u ) - log(l + e v ~ u ), 
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u x =e Ul ~ u + e Ul - v , 
Uy = e «-i-« + e «-i-« 

a = l + e v - u , H-- 
u x = e Ul ' v + e u ~ v , 



Vx = e «-«-i +e u - v -\ 



Ml— U 



i? 



«— ui 



u„ = e *+e" v y = e ul " + e 



U— U_l _|_ i>— « 

a = 1, F = -e" 



II 



II „ 



U-\—U 



V\—U _|_ v — u 

# = e «i-u + e «- Uj 

U— V—l I it— D 



3 «— «i 



g«_l— U _ 2 y git — _ J ' 

a = _ e ^« H = e Ul ~ v + ^e 2{u ~ v \ R = log(l - e"" 1 ' 1 ). 

7.2 Hyperbolic PDE systems 

o^xj/ = o, u u x u y + 2a(u x — u y ) + aa„ — a„a, 



■<\v 



a v v x v y - 2a(v x - v y ) + aa u - a u a, 



4a 



TO 

(n) 

TO 

TO 
WO 



(itj, - l)u x 
Uxy = — h ltj/(« - U), 



m 



m — M 
Ita-^ - 1) 



■ni 



^ y ' +u y (l + e u - v ), 
l + e u-v v\ » 



u 



a 



1 + e u ~ v 



+ l + e l 



it 



X y U X Uy -j- 26 Uy 1, 

7.3 Ruijsenaars— Toda lattices 

u xx = (u 2 x + r(u)) 
Uyy = {u 2 y + r{u)) 



J xy 



J xy 



(Vy ~ 1)V X 

u — V 

V X (Vy ~ 1) 

1 + e u ~ v 



- Vy(u ~ V), 



1 + e u ~ v 



l-e l 



.'/- 



v x v y - 2e u v v y + 1. 



Ul, x 


U-l,x 


T(a) 


a 


«l, v 




T(a) 


a 



+ i^log(r(o)a)) -^(u), 
+ i^bg(T(a)a)) -ir'(u), 

(n^-Uy - r - u_i) + (it x + u y )(a + -a w _ 1 (ui - tt_i)) = 
k a = a(u,it_i), 4r(u) = 2aa u _ lU _ 1 - a^, 

mi + 2m — u_ i I , 



u 



Ml — U U — U-\ 



yy 



-Uy(Uy ~ 1) 



Ml — U U — M_i 

k U x Uy = (1 - Uy){u\ ~ U)(U ~ M_l), 

Ul, x U-i x 



Ul—U _|_ gM-W-i 



l + e u- U1 l_|_gU_l-U 



M 



yy 



[ = K - l)(e" + l)^"-"- 1 + 1), 



(H 2 ) 

m 



(i?2 



(i?3 
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u xx = u x {m )X - u- 1>x - e Ul u + e u u ~ 1 ), 

Uyy = Uy{U- lt y - U^y - e"" 1 "" + e""" 1 ), 



( U 1, X U X 



U X U-1 



u 

, U X Uy 



l + e u- Ul ]_ + e u_i-u' 



l + e u- Ul i + e «_i-u' 
= (e" 1 - n + l)(e u -"- 1 + l), 

e* 



2(ui— u) _|_ g 2(u— 



"1/2/ ~~ J/ 
. U x - l/u. 



(i? 4 ) 



(R 5 



u _|_ gU— u— 1 



Acknowledgements 

This research was supported by the RFBR grant # 04-01-00403. 



[1] Ruijsenaars S.N.M., Relativistic Toda system, Preprint Stichting Centre for Mathematics and Computer 
Sciences, Amsterdam, 1986. 

[2] Ruijsenaars S.N.M., Relativistic Toda systems, Comm. Math. Phys., 1990, V.133, 217-247. 

[3] Bruschi M., Ragnisco O., On a new integrable Hamiltonian system with nearest-neighbour interaction, 
Inverse Problems, 1989, V.5, 983-998. 

[4] Suris Yu.B., Discrete time generalized Toda lattices: complete integrability and relation with relativistic 
Toda lattices, Phys. Lett. A, 1990, V.145, 113-119. 

[5] Adler V.E., Shabat A.B., On the one class of the Toda chains, Theor. Math. Phys., 1997, V.lll, 323-334. 

[6] Adler V.E., Shabat A.B., Generalized Legendre transformations, Theor. Math. Phys., 1997, V.112, 935-948. 

[7] Adler V.E., Shabat A.B., First integrals of the generalized Toda lattices, Theor. Math. Phys., 1998, V.115, 
349-358. 

[8] Marikhin V.G., Shabat A.B., Integrable lattices, Theor. Math. Phys., 1999, V.118, 217-228. 

[9] Adler V.E., Marikhin V.G., Shabat A.B., Canonical Backhmd transformations and Lagrangian chains, 
Theor. Math. Phys., 2001, V.129, 163-183. 

[10] Ragnisco O., Santini P.M., A unified algebraic approach to integral and discrete evolution equations, Inverse 
Problems, 1990, V.6, 441-452. 

[11] Shabat A.B., Yamilov R.I., Symmetries of nonlinear lattices, Leningrad Math. J., 1991, V.2, 377-400. 

[12] Adler V.E., Yamilov R.I., Explicit auto-transformations of integrable chains, J. Phys. A: Math. Gen., 1994, 
V.27, 477-492. 

[13] Yamilov R.I., Symmetry approach to the classification from the point of view of the differential-difference 
equations. Theory of transformations, Doctoral Thesis, Ufa, 2000. 

[14] Adler V.E., Shabat A.B., Yamilov R.I., Symmetry approach to the integrability problem, Theor. Math. 
Phys., 2000, V.125, 1603-1661. 

[15] Adler V.E., Discretizations of the Landau-Lifshitz equation, Theor. Math. Phys., 2000, V.124, 897-908. 

[16] Hirota R., Nonlinear partial difference equations. II. Discrete-time Toda equation, J. Phys. Soc. Japan, 
1977, V.43, 2074-2078. 

[17] Suris Yu.B., Bi-Hamiltonian structure of the qd algorithm and new discretizations of the Toda lattice, Phys. 
Lett. A, 1995, V.206, 153-161. 

[18] Adler V.E., On the structure of the Backlund transformations for the relativistic lattices, J. Nonlinear Math. 
Phys., 2000, V.7, 34-56, [nrImSI/0001072| 

[19] Adler V.E., Suris Yu.B., Q4: Integrable master equation related to an elliptic curve, Internat. Math. Res. 
Not, 2004, V.47, 2523-2553, [nlin.SI/0309030| 

[20] Zhiber A.V., Ibragimov N.H., Shabat A.B., Liouville type equations, DAN SSSR, 1979, V.249, 26-29. 

[21] Sokolov V.V., Zhiber A.V., On the Darboux integrable hyperbolic equations, Phys. Lett. A, 1995, V.208, 
303-308. 



On the One Class of Hyperbolic Systems 



17 



[22] Zhiber A.V., Sokolov V.V., Exactly solvable hyperbolic equations of the Liouville type, Uspekhi Mat. Nauk, 
2001, V.56, N 1, 63-106. 

[23] Yamilov R.I., On classification of discrete evolution equations, Uspekhi Mat. Nauk, 1983, V.38, N 6, 155-156. 

[24] Adler V.E., Shabat A.B., Dressing chain for the acoustic spectral problem, Theor. Math. Phys., 2006, V.149, 
1324-1337, |nlin.SI/0604008| 

[25] Nijhoff F., Hone A., Joshi N., On a Schwarzian PDE associated with the KdV hierarchy, Phys. Lett. A, 
2000, V.267, 147-156, [solv-int/9909026j 

[26] Svinolupov S.I., Yamilov R.I., The multi-field Schrodinger lattices, Phys. Lett. A, 1991, V.160, 548-552. 

[27] Merola I., Ragnisco O., Tu G.-Z., A novel hierarchy of integrable lattices, Inverse Problems, 1994, V.10, 



1315-1334, solv-int/9401005 



